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Abstract 



Stimulated by the results of NMR experiments with superfluid 
%^J 3 He in "nematically ordered" aerogel pQ we report on the results of 

^^ . phenomenological analysis of stability of different phases of superfluid 

CS \ 3 He subjected to a strong homogeneous uniaxial anisotropy. On a 

basis of this analysis we suggest a form of the order parameter for the 

new ESP2 phase observed in the quoted experiments. In the weak 
^^ [ coupling limit the suggested order parameter approaches that of the 

axi-planar phase. We discuss a possible experimental check of the 

proposed identification of the new phase. 

1. Aerogels opened a possibility to impose a global orbital anisotropy 
on the p-wave Fermi superfluid - 3 He [2] . In the mean field approxima- 
tion deformed aerogels couple to the orbital part of the order parameter in a 
way analogous to Zeeman coupling of the spin part to magnetic field. Corre- 
sponding term in the density of free energy is proportional to KjiA^A*^ . The 
complex 3x3 matrix A^ is the order parameter of 3 He, averaged over dis- 
tances larger than the distance between the strands of aerogel £ a . The index 
jjl stands for spin projections and j - for orbital. Real symmetric traceless 
tensor Kji describes effect of the anisotropy. If anisotropy is uniaxial and 



z is the symmetry axis Kji can be written as k xx = K yy = — k , k zz = Ik . 
The uniaxial anisotropy splits the transition temperature for orbital triplet 
{1=1) in two different temperatures- for l z = and for l z = ±1, favoring 
one or the other state, depending on the sign of k . Uniaxial compression 
( k > ) favors l z = ±1 projections. It increases stability of the ABM phase 
and orients its order parameter. Physical manifestations of this orientation 
were discussed in the literature [3J H]. Uniaxial stretch ( k < 0) favors l z = 
projection, which taken alone can not form the orbital part of the ABM order 
parameter. Aoyama and Ikeda [5] predicted that this type of anisotropy will 
create just below T c a region of stability of the polar phase. The width of 
this region depends on the absolute value of k , i.e. on a rate of anisotropy 
of aerogel. Possibility of stretching of silica aerogels is limited. To obtain an 
appreciable anisotropy in the experiments [T] a new type of aerogel was used 
[6], [7]. Its strands are nearly parallel to one direction (anisotropy axis), so 
that it was dubbed as "nematically ordered". Anticipating we remark that 
even for that new aerogel the estimated from the experiment value of k is of 
the order of 10 ~ 2 and the anisotropy has strong effect on the phase diagram 
of 3 He in aerogel only near to its transition temperature T ca . 

The authors of the paper [1] investigated superfluid phases of 3 He in 
aerogel by CW and pulsed NMR methods for different orientations of the 
DC magnetic field with respect to the anisotropy axis (taken as z-axis 
in what follows). Several phases with different NMR-signatures were ob- 
served and the phase diagram of 3 He in " nematically ordered" aerogel in the 
temperature-pressure coordinates was constructed. The results, obtained on 
cooling from the normal phase confirm the scenario, suggested by Aoyama 
and Ikeda. Helium first enters the Equal Spin Pairing (ESP) state which 
can be interpreted as following the root from the polar phase via the second 
order phase transition to the ABM-phase with polar distortion. The latter 
on further cooling jumps to the low-temperature phase (LTP). This phase 
has the NMR properties of the BW phase with the polar distortion. 

Surprising results are obtained on warming from the LTP. On its ap- 
proach to the T c helium jumps to the ESP-phase which has different NMR 
frequency shift from one observed on cooling. The new phase was denoted as 
ESP2 and it is not properly identified. To clarify this situation we carried out 
an analysis of possible phase transitions in uniaxially stretched aerogel phe- 
nomenologically, within the Landau theory of phase transitions. Superfluid 
3 He in aerogel can be considered as a uniform media if the Ginzburg and Lan- 
dau coherence length £(T) is greater than the average distance £ a between 
the strands of aerogel. In terms of the zero temperature coherence length £ 
this condition can be reformulated as a restriction of the temperature interval 
where the state of helium can be characterized by average order parameter: 



(T - T ca )/T ca < (^o/^a) 2 • For the low pressure data of Ref.JI] (P < 6.5 bar) 
the Landau theory can be applied down to w .85T ca . This region includes 
the most interesting phase transitions. For the high pressure data (P > 12 
bar) £o <C £ a and the interval of applicability of uniform approximation is 
several times smaller. Outside of this interval it can be used only as a quali- 
tative guidance because the state becomes essentially nonuniform and local 
fluctuations of the order parameter can render significant contribution. 

2. Within the mean field approximation the average thermodynamic 
potential of the p-wave superfluid in a globally anisotropic aerogel has the 
following form: 

$ s = $ n + JVe//[(T<J;i + KjfiAnA^ + -{^A^A^Al^ + f3 2 A M A* M A vl Al l + 

PzA^jAvjA^A^ + p^A il jA V jA v [A i)X + p^A^jA^A^A^)] (1) 

Here t = (T — T c0 )/T c0 is the dimensionless temperature counted from the 
transition temperature at k — 0, N e ff denotes an "effective density of 
states". The values of the phenomeno logical coefficients /3 1; .../3 5 depend on 
pressure and properties of aerogel. The values of @i, ...p 1 *, derived from the 
BCS theory are referred as the weak coupling limit. When normalized to fa 
these values are /?i, /3 2 , /3 3 , /3 4 , /3 5 = /3 2 (-1/2,1,1,1,-1) [$]. Experimental data 
indicate that real values of these ratios even for the bulk 3 He, may deviate 
from their weak coupling values for 20 % - 30 % [5]. The deviations are 
important at determining regions of stability of different phases. 

For uniaxially anisotropic aerogel the second order terms in Eq. (1) have 
the form (T+2K)A flz A* flz + (T-K)(A flx A*^ x + A liy A^ y ) . If « < the superfluid 
transition takes place at r = —2k . Below this r the order parameter is a 
complex spin vector A^ z = AqA^J with A® normalized by the condition 
AJJ'JAjJ']* = 1 . Minimization of the free energy with that order parameter 
renders 

A 2 t + 2k (2 ) 

Pl5\-n-Hz Apz \ + P234 

Here the conventional shorthand notations /3 15 = /3i + /3 5 etc. are used. 
All data indicate that in 3 He /3 15 < , then the maximum gain of energy 
is reached at |A( JA( j| 2 = 1. Together with the normalization condition 
it renders A® = exp(iip)d^, where d^ is a real spin vector. With the 
account of the orbital part it reproduces the order parameter of the polar 
phase A°j = Aoexp (if) d^rrij , where rrij is a unit vector in z -direction in 
agreement with the Ref.[Sj. For the overall amplitude we have from Eq. (2): 
A2 = -(r + 2 K )//3 12345 . 



3. To find further possible phase transitions we represent the order pa- 
rameter as A^j = A ^ +a^j , where a^ is a small increment, and expand the 
change of the thermodynamic potential $ = $ s — $„ in powers of a^j . All 
experimentally observed transitions take place at \r\ < .15, where Eq.(l) is 
still a good approximation for $ . Temperature of the transition is deter- 
mined by the second order terms: 

*(A ltj ,A; ij ) = *(A%,A%)+ 

i \ d 2 § <9 2 $ , d 2 § ,1 

2 J dA M dA vl "*"« + dA^dAU a "* a « + dA^dAU a * a « / " {6) 
At the transition the linear equation 

<9 2 $ d 2 § 

dA-M;^ 1 + di-dAj^ ~ ° (4) 

together with its complex conjugated equation acquires a nontrivial solu- 
tion, which breaks the symmetry of the polar phase. To select essential 
solutions we have to impose the orthogonality condition A°*-a w - = , or ex- 
plicitly: d^rrija^j = . Because of the degeneracy of A Q ■ with respect to 
orientation of the spin part d^ even more restrictive condition can be im- 
posed vdja^j = . Using the freedom in determining the overall gauge we 
require that for the polar phase exp(i<p) = 1 , then A • is real. Eq. (4) can 
be rewritten as two equations separately for the real and imaginary parts 
of a^j = a R j + ia 1 ^ . Both equations have the form K '^a^ = with 

Knvji — A^SfivSji + B^S^rrijmi + C^d^Sji + D^d^d^rrijini . Because 
of the imposed conditions the terms with B R)I and D R)I do not contribute 
to the equations for a^ . The resulting equations have the form: 

(A R > T V + C^d^a^f = 0. (5) 

Let dft, e M , f^ be orthogonal basis in spin space. Taking projections of Eq. 
(5) on each of these vectors we arrive at the following possibilities: 1) A R)I + 
C^ 1 = , d^f ^ or 2) A R > T = , e^a*f ^ , and f^a*f ^ . In both 
cases the orbital part is orthogonal to rrij . Coefficients A R)I and C R)I can 
be expressed in terms of r, k and fli, ...{$$ with the aid of Eqns. (1),(2). 

First consider the imaginary part. The sum A 1 + C 1 = r — k + (/?245 — 
/3i2)Aq turns to zero at t — Ta, ta = t(3/?245 — fiiz)/1fiiz . The emerging 
component is 

a 1 ^ = iAidpUjt (6) 

where rijl is a unit vector, orthogonal to rrij and Ai is a real amplitude. 
Together with A ■ it forms the order parameter of the ABM-phase with a 



polar distortion A w = d^(A mj + AiUji) , which was discussed in the Refs. 
[5J, pp. The temperature dependence of A and A 1 at r < t a is found 
by minimization of $ s (Eq. (1)): A? = -g£, Ag = __|_ _ igi . T n 
comparison with the polar phase the new phase has lower thermodynamic 

potential. The gain is &a — $ p = — 2 „ ^" (r — r^) 2 . The condition 2) for 

imaginary part A 1 = is not satisfied at small r . 

For the real part the sum A R + C R = —3k is finite at all r . It means 
that the real part of the projection of a MJ - on d M is absent. On the contrary 
A R = r — k + /3i 2 Aq turns to zero at r = r B , r B = «(1 + 3/3i 2 //3 3 4 5 ) . 
Below this temperature finite projection of a R , on e^ and / M can exist. 
With the aid of rotations around d^ and m^ the 2x2 real tensor in a space 
orthogonal to these vectors can be transformed to the form 

a*=A 2 e li l j + A 3 f fl n j , (7) 

rij and lj are two mutually orthogonal vectors, forming together with rrij 
basis in the orbital space, A 2 and A 3 are real amplitudes. Minimization of 
$ s at r < t b renders: A 2 = A 3 , then A - + a R j forms the order parameter 
of the BW-phase with a polar distortion, also discussed in Refs. [5], [I]. 
Temperature dependence of the amplitudes is given by A 2 , = — 3l3 T ~+a , 

A o = -fc + A 2 ^th the energy gain $ B - % = -^-Jm^—^t - r B f . 

On cooling from the polar phase the perturbation with the higher tran- 
sition temperature r occurs first. According to the experiment [T] this is 
the ESP-phase. Of the two considered possibilities only a*j meets this 
requirement. It means that ta > t b and the ABM-phase with a po- 
lar distortion develops below ta ■ Repetition of the same argument with 
ApJ = Aod^Tiij + iAidfjij shows that no further continuous transition oc- 
curs on cooling within the limits of applicability of the expansion Eq. (1). 
The more advantageous at low temperatures B-phase with a polar distortion 
occurs via the first order transition in agreement with the experiment [1] and 
with the theoretical analysis of Ref. [5]. Since the sequence of the transitions 
depends on the sign of the difference ta~t b = 3k/3i 2 345(1/2/3i 3 — l//3 345 ) we 
conclude that for 3 He in nematically ordered aerogel /3 345 < 2/3 i3 . 

4. The case t a = t b = r AB is special. It is realized in particular for the 
weak coupling values of the coefficients /?i, ...f3 5 . In this case the emerging 
increment of the order parameter can be searched as a linear combination 
of a • and a R , so that the order parameter just below the transition has a 
form: 

A^j = Aod^rrij + iAid^rij sin 9 + lj cos 6) + A 2 e fl l j + A^f^nj, (8) 

where 6 is the angle between nfj and lj . Minimization of thermodynamic 
potential (1) with the weak coupling values of (3i, .../3 5 renders both states 



considered above plus one more ESP-state: 

A* p = Aod^rrij + iA^Uj + A 2 e fi l j , (9) 

which is a particular form of the axi-planar phase [10j,[llJ. For the am- 
plitudes A ,Ai,A 2 we have: Aj = ^ (^ - r) , A 2 2 = 3A? , Ag = 

-i^( 2r + 7k ) and the energy gain ($ p - $ AP )/N eff = ^(r - r A f - i.e. 
exactly the same as that for the distorted ABM-phase. At ta = tb both 
ESP phases are less advantageous than the distorted B-phase, so that the 
B-like distortion develops at r < tab ■ 

Deviation of the coefficients /3 S from their weak coupling values results 
in splitting of the transition. If the splitting is small i.e. (ta — Tb) <C ta , 
the axi-planar order parameter Eq. (9) with possibly corrections of the order 
of {ta — t~ b )/ta in a region \r — ta\ ~^> (ta ~ t b) renders an asymptotic 
form for a local minimum of the thermodynamic potential, which can be 
reached via the first order transition or a sequence of transitions. To find out 
properties of the new phase we have to substitute the order parameter (8) 
in the expression for the thermodynamic potential Eq. (1) and minimize it 
over 5 free parameters: 9, Ao, Ai, A 2 , A3 . Extremum over 9 is reached at: 
1) A\ = , 2) Ajj = A§ , 3) sin6> = or cos6> = (these two possibilities 
are equivalent up to notations). The possibilities 1), 2) correspond to the 
distorted BW-phase. The most interesting for the present discussion is the 
third possibility. For definiteness we take cos 8 = 0. If all A;? are finite they 
have to obey the following system of linear equations: 

/3i234 5 A2 + (/3 245 - Pn)A\ + (3 12 A 2 2 + (3 l2 A\ = -(r + 2k), (10) 

(/3 245 - /3ia)Ag + /3i 2 345A? + (f3 2 - 13,) Aj + (/3 234 - /3is)A§ = -(r - «), (11) 

(3 l2 Al + (/3 2 - p x )A\ + /WsAo 2 + (5 l2 A\ = -(r - «), (12) 

P12AI + (/3 234 - /?i 5 )A? + {3 12 A\ + /3 12345 A2 = -(r - k). (13) 

The lack of precise knowledge of the values of /3 -coefficients introduces ambi- 
guity in analysis of these equations. Here we remark only on some properties, 
which follow from the assumption that the strong coupling corrections are 
small. Combination of Eqns. (12) and (13) renders the relation: 

((3 M -(3 5 )Al = p M5 (AJ-Al). (14) 

As a result if A^ 7^ A 2 > A| . In a vicinity of the singular point ta = t~b of 
particular interest are combinations of the coefficients /3 , which tend to zero 
in the weak coupling limit. The expression for the splitting ta — t~b i n the 



leading order depends on two such combinations /3 45 = e/3 2 and 2p\ + (3 3 = 



V$2 



\JA ~ T B 



K£ — V) 



Combination of Eqns. (11) and (12) renders 



[&4 - (2ft + /3 5 )]A 2 = [{u-e) A 2 + (z/ 



e)(Al 



A?)]A. (15) 



At finite e and z/ A| 7^ so that the new phase is strictly speaking non-ESP 
and not axi-planar. In a principal order on e<l and i/<1: 



A 2 ; 



1 



16ft 



(3i/ - e)r + (9i/ 



. K 



(16) 



Presence of the term, proportional to r in this expression means that the 
order parameter of the new phase remains distorted even in the limit |r| ^> 
|«| . The distortion is relatively small only because e and v are small. 

5. Now we compare the expected NMR properties of the new phase 
with the properties of the ESP2-phase of Ref. pQ. In these experiments 
different phases are characterized by their transverse NMR frequency shifts 
as functions of the tipping angle j3 and the angle \i between the direction of 
the d.c. magnetic field and the anisotropy axis z The calculated dependence 
of the shift on these angles for the order parameter (8) with cos 9 = in the 
two-dimensional Larkin-Imry-Ma state is given by the formula: 



XCo^Aw = const. (2|A | 2 — |Ai| 



IA I 2 



sin 2 /1 



(7cos/3 + 1) — cos/3 



sinV 
cosp — (5 cosp 



+ |A 3 | 2 cos/3 2-3sinV 



+ 



(17) 



Here the proportionality coefficient is the same for all phases and it is deter- 
mined by the properties of the normal phase, ojl is the Larmor frequency 
and x ~ t ne maximum principal value of a tensor of the magnetic suscepti- 
bility of the considered phase. As a zero order approximation in e and v we 
use the order parameter A^ with temperature dependencies of A , Ai , 
A 2 as given in the paragraph, following Eq. (9). Then 



X^lAwap 



1 



19 \ sin 2 /1 

2v r+ rv cos/3+ ^6- 



K 



(7cos/3 - 5)r + -(205 cos (3 - 23) 



For comparison, in the distorted ABM-phase (ESPl-phase of Ref.pQ): 



X^lAwabm 



1 / 19 

2 T + Y K 



cos/3 



(sin 2 jx) (5 cos (3 — 1) 



(19) 



When magnetic field is parallel to z ( fi =0) the calculated transverse NMR 
shift for the axi-planar phase coincides with that for the ABM-phase. In the 



experiments [I] the shifts for ESP1 and ESP2 phases do not coincide, but the 
difference is not big and it can be ascribed to the corrections of the order of 
e and v . More definitive experimental check of the proposed identification 
would be a measurement of the NMR shifts at the perpendicular orientation 
of the field to the anisotropy axis ( sin 2 \i =1). For this orientation two shifts 
have to be different - in the distorted ABM-phase CW shift is zero, while in 
the axi-planar phase 

XU L Au A p -(t--kJ. (20) 

The ratio of the shift for sin 2 // =1 to that for sin 2 /i =0 in this phase has to 
be J 2 t+i9k) a ^ (^ r — 5k) < 0. The data for the ESP2 phase at sin 2 /i=l 
are not yet available. 

In conclusion, we suggest to identify the ESP2 phase, observed in the 
experiments [1] as the phase with the order parameter Eq. (8) at cos 9 = 0. 
In the principal order on small deviations of coefficients (3 S from their weak 
coupling values it reduces to the axi-planar form A^ given by Eq. (9). The 
proposed identification can be checked in NMR experiments. The performed 
analysis of possible phase transitions in the Ginzburg and Landau region can 
be used also for extraction of values of coefficients /3 S from the experimental 
data with a subsequent quantitative interpretation of the data. A more 
extensive experimental investigation in particular in the low pressure region 
would be useful for realization of this program. 
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